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Abstract 

The cross section for the reaction e~e~ — > i~t~ [t = n,r) is calculated in models with 
heavy Majorana neutrinos mediating lepton number violating amplitudes at the loop level. The 



contributing four-point functions are evaluated exactly (numerically) taking into account the full 

CM ' 

propagator dependence on external momenta, thereby extending to the energy range of interest 
for the next linear colliders an earlier approximate low energy calculation. The amplitude shows a 



non-decoupling behaviour relative to the heavy Majorana neutrino masses, but due to the stringent 
bounds on heavy-light mixing the signal cross section attains observable values only for the less 



o 

cm 

constrained r signal. The cross section induced by lepton number violation in the SU(2)i doublet 



sneutrino sector of supersymmetric extensions of the standard model is constrained by the upper 



limits on neutrino masses and probably too tiny to be observable. 
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The process e~e~ — > £~£~ (£ = /x, r), which violates the L e and L £ lepton numbers, 
is forbidden in the Standard Model (SM) due to exact lepton number conservation to all 
orders of perturbation theory. Its observation at a next generation linear collider with 
center of mass energy = 500, 800, 1000 GeV may be possible only if there is new 
physics that can trigger it. The signature is clear and practically free from SM background. 
In literature it was studied: (i) in the context of models with gauge bileptons jlj], where 
the final state is reached through tree level s-channel annihilation into a gauge bilepton and 
subsequent decay, (ii) in the context of mixing models where the reaction proceeds through 
a loop (box diagram) with heavy Majorana neutrinos and W~ gauge bosons as virtual 
particles running in the loop |2j, and (Hi) in supersymmetric scenarios where sneutrinos and 
charginos instead of neutrinos and charged bosons are exchanged {jj . The aim of this paper 
is: (i) to improve and extend the calculation of Ref. Q (which was essentially a low energy 
calculation) in order to provide predictions in the energy range of interest for the next linear 
collider project, motivated by the observation that all diagrams of the process, see Fig. ^a- 
d), have a threshold singularity at y/s = 2Mw where the amplitude develops an imaginary 
part givinga boost to its absolute value (see the well known example of photon-photon 
scattering |4J), taking also into account experimental bounds on effective mixing angles not 
considered in [2J|; (ii) to make a realistic calculation for the sneutrino case: the cross section 

n 

was estimated so far assuming eV scale sneutrinos [3J. 

We assume that heavy Majorana neutrinos (mass eigenstates) couple to the standard 
model charged currents through heavy-light neutrino mixing. This is the simplest way to 
obtain lepton number violating processes. Consider the box diagrams depicted in Fig. ^a- 
d). As in we use the 't Hooft-Feynman gauge: there are graphs with WW, 00 and <fiW 
exchange, <fi being the Goldstone boson. The lagrangian of interest is l^ : 



i,Ni v2 



M Ni - l + 7 5 rr , , 
-ip e — - — U eNi ijj Ni (j) 



+ h.c. (1) 



M w 2 

where £ = e, /x, r and are the elements of the mixing matrix of the heavy mass eigen- 
states M^- labeled by the index N%. Neglecting the masses of the external particles allows 
to simplify the calculation of the amplitudes that can be expressed in terms of: (i) the 
Mandelstam variables s, t and u; (ii) the spinor products of light-like momenta (see Ref. ^| 
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FIG. 1: In (a-d) the Feynman diagrams, in the 't Hooft-Feynman gauge, contributing to e~e _ — > 
i~ i~ [l = fj,, t) via heavy Majorana neutrinos. In (e) the corresponding SUSY diagram is given 
which arises in models with lepton number violation in the sneutrino sector. In (f) the choice of 
the running momentum is given. The momenta ki (i = 1, 2, 3) for the decomposition of the tensor 
integrals within the LOOPTOOLS notation are: k\ = pi, k 2 = Pi + P2, ^3 = Pi + P2 — P3 = Pi- 

and references therein) 



S(Pa,Pb) = U + {p a ) • W_Ofc) 
T(p a ,Pb) = «-(Pa) ■ U + (p b ). 



(2) 



which obey the relation: \S(p2,pi)T(p 3 ,p4)\ 2 = s 2 ; and (Hi) the scalar D and the tensor 
rank-1 (D^) and rank-2 (D^) four-point functions jf|. The corresponding amplitudes are 
found to be: 

4 i 

' E (U: N U, N] fM N M Nj 

Ni,Nj 

(3) 



M 



WW 



JL 



1 



\V2J (4tt) NifN . 
x AS(p 2 , Pl )T(p3, p 4 ) [D (s, t) + D (s, u)} 



Mr = A= 



1 Vfrr rr V M%M* Nj 



M W M W 



x M Ni M N] S(p 2 , Pl )T(p 3 ,p 4 ) [D Q (s,t) + D (s,u)] 



(4) 
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x S(p 2 ,pi)T{p 3 ,p 4 ) [4(D 00 (s,t) + D 00 (s,u)) 

- 2(tG(s,t)+uG{s,u)) 

- 2(tV(s,t)+uV(s,u))} , (5) 



M 4>W 19 1 1 V /m tj , 2 M Nl M Nj 



x S(p 2 , Pl )T(p 3 ,p 4 ) [A(D 00 (s,t) + D 00 (s,u)) 

- 2(tG(s,t)+uG(s,u))] . (6) 

The numerical computation of the four-point functions was performed using the LOOPTOOLS 
[z] software where the following notation is used as regards the expansion of the rank-1 and 
rank-2 tensor functions: = £? =1 A (h)^ D^ = g^D 00 + Yhj=i D ij (^)m (^i)xo ^ bein § 
sums of external momenta running in the loop as explained in Fig.^f). Within this notation 
the form factors G and V appearing in Eqs. ©-© are given by 

G = D 22 + D 23 + D 12 + D 13 , 
V = 2D 2 + D X + D 3 + D . 

Terms depending both on (s,t) and (s,u) appear because the identical fermions in the final 
state require proper anti-symmetrization of the amplitudes. Defining xw = sin 2 9w and 
Xij = M%. ./Myy, the differential cross section is easily found to be: 

da 1 / a 



(—) {Kis^u^s, (7) 



dcos9 2567T \xw / 
where K is given by: 

k = E (u: Ni u^) 2 v^]{m 2 w (1 + x -f) 

x [D (s, t) + D (s, u)] + 2 (D 00 (s, t) + D 00 (s, u)) 

us-it ,\ ^/ M \tV(s,t) +UV(S,U)]} /nN 

- [tG(s,t)+uG(s,u)\- 1 ' 2 K - 2 - J1 \. (8) 

To obtain the total signal cross section a to t, Eq. is integrated numerically over the 
scattering angle in the center of mass frame. As stated above, similar formulas were obtained 
in P| using the approximation where all external momenta in the loops are neglected relative 
to the heavy masses of the gauge bosons and Majorana neutrinos, enabling to carry out the 



loop integration analytically. The formulas thus obtained are well known in the literature [8( 
and the final cross section, which depends only on scy and the mixing coefficients, grows 
linearly with the center of mass energy squared, s. This approximation for the four-point 
functions is good at low energies, such as in decay processes of heavy mesons, or when 
\/s << M, M being the highest mass running in the loop. In addition the linear growth 
with s would break unitarity, therefore in order to make quantitative predictions with the 
correct high energy behaviour, the four-point functions full dependence on the external 
momenta has to be considered. Theoretically, according to the 'Cutkosky rule', one expects 
an enhancement at ~ 161 GeV ~ 2Mw, the threshold for on-shell WW gauge boson 
production, at which the four-point functions develop an imaginary part. In Fig. Efa) the 
ratio R a = Gtot/&o of the integrated total cross section a to t to <7o, the cross section of the 
low energy calculation of Ref. 2j, is plotted for sample values of the Majorana masses. The 
enhancement due to the threshold singularity of the loop amplitude is more pronounced for 
values of Majorana masses close to Mw and is drastically reduced increasing ~ to 
O(TeV). As R a ^ 1 as yfs « My/ in all the cases, the agreement of our full calculation 
with the result of Ref. |2| in the regime of low energies is evident. 1 

The threshold effect appears to be quite spectacular only for values of Majorana masses 
that correspond to cross sections too small to be measured even at a next linear collider. 
In Fig. 12(b) the effect of the threshold singularity in the loop integral is shown reporting 
absolute cross sections for a particular choice of Majorana masses: = 150 GeV and 

MjVj = 450 GeV. The low energy approximation (dashed line) obtained neglecting external 
momenta in the loop is inadequate when the energy of the reaction increases to values 
comparable with the masses. Increasing the energy, after reaching a maximum, the cross 
section starts to decrease until the asymptotic behaviour 0(1/ s 2 ) of the loop integral K is 
reached. This happens for every value of heavy Majorana neutrino masses and we checked 
numerically that, as expected, for higher masses the asymptotic regime is reached at higher 
values of yfs. In fact from Fig. El we note that the cross section grows with increasing HMN 
masses. The main contribution comes from the graph with two Goldstone bosons since their 



1 It should be mentioned that the agreement, at low energies, of our Eq.(7) with Eq.(9) of Ref. 2] is up to a 
factor of 4. We have contacted the author of Ref.|2l on this matter and he agrees with our Eq.(7). That 
is, Eq.9 of Ref. [3] should be divided by 4 (in Refjjl] the average over initial spins was left out 0) and 
then for energies y/s << Mw it coincides exactly with our Eq.(7). 
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(a) (b) cr (fb) 




Vs (GeV) 

FIG. 2: In (a) the ratio R a is plotted as a function of ^fs the energy in the center of mass 
system: solid line, M N . = M Nj = 100 GeV; short-dashed line, M N . = 150 GeV, M Nj = 450 GeV; 
long-dashed line, M Ni = M Nj = 500 GeV; dot-dashed line, M Ni = M Nj = 1 TeV; dotted line, 
Mtv ; = iWjVj = 3 TeV. In (b) the absolute values of the cross sections are given for a particular 
choice of Majorana masses = 150 GeV, Mj^. = 450 as function of the energy in the center of 
mass frame. The solid line is obtained integrating Eq. (J7J) while the dashed line is based on Eq. (9) 

of Q. 

coupling is proportional to M^. Moreover the chiral structure of the coupling selects the 
mass term in the numerator of the Majorana neutrino propagators. When these masses 
are much larger then the other quantities, the amplitude scales like M^Mff./M^M^. ~ 
M/v-M^ , i.e. is proportional to the square of the heavy masses. This fact is the well known 
non decoupling of heavy fermions in theories with spontaneous symmetry breaking (similarly 
in the SM the top quark gives sizable radiative corrections owing to its large mass and a 
quadratic non decoupling). 

Heavy Majorana neutrinos naturally appear in extensions of the SM with right-handed 
neutrinos which generate light neutrino masses through the see-saw mechanism. The scale 
of the masses Mr is of order 10 9-12 GeV with very small heavy-light mixing, L^j~M^j, 
and the cross section will be suppressed by inverse powers of these masses, recovering so the 
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FIG. 3: Total cross sections as function of the center of mass energy, y/s. The value of the mixing 
coefficients are discussed in the text. In part (a) the solid curve referes to the case of e~e~ — > t~t~ 
with = M/Vj = 3 TeV, while the dashed line referes to (e~e~ — > /x/i) with the same values of 
Majorana masses. In (b) the Majorana masses are changed to somewhat lower values: = 1 



TeV, Mjv, = 3 TeV. 



decoupling limit that is natural in the see-saw framework. 

An interesting scenario is the model recently proposed in Ref. [24], where an attempt 
is made to construct light neutrino masses with a see-saw mechanism and no new physics 
beyond the TeV scale. This is achieved by adding a new Higgs doublet, relative to the 
SM, whose neutral component develops a naturally small vacuum expectation value, u ~ 1 
MeV, so that m v = m 2 D /M N = f 2 u 2 /M N ~ 1 eV if M N ~ 1 TeV and / ~ 1, with / being 



lich is too small 
that a charged 



the Yukawa coupling. But the heavy-light mixing is ~ fu/M^ ~ 10 -6 , w 
to have phenomenological consequences. Further it was shown in Ref. [2£ 
Higgs boson of this model must be heavier than 50 TeV, in order to satisfy the experimental 
bound of the /i — >• 3e decay. So this model does not comply with the constraints from 
non-observation of lepton flavour violation. 

More interesting from the phenomenological point of view is the case in which HMN 
have masses in the TeV range with non negligible mixing. Mass matrices that satisfy this 
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condition can be built, using experimental constraints on heavy-light mixing (including those 
from /3/3o v This is achieved imposing relations among the elements of the neutrino 

mass matrix in a way that the mixing is decoupled from mass relations and is bounded 
only by data Q.Q- Independence of the mixing matrix from the mass Nation and the 
consequent possibility of violating the Appelquist-Carazzone theorem |14j have led many 

0,0,0, 



aut 
16, 



rors to study HMN contributions to rare processes like 



17 



s 



e7, fi 



e 1 e e 



[rsl ll^ . However it was recently argued in Ref . [11| that even if such a situation is 
not still ruled out by present data on neutrino oscillations, it requires extreme fine tuning 
among the elements of the Dirac mass matrix mo and those of Mr. Keeping this in mind, 
we can nonetheless explore the phenomenological consequences of such a scenario. As was 



done in Ref. 
mixing [ill IsC 



we take the following experimental upper bounds on effective heavy-light 
211= 



]T \U eNi \ 2 < 0.0027, 

Ni 



(9) 



< = £|£W 2 < 0.005, 

Ni 

< = El^l 2 <°- 016 > 

Ni 

and allow the heavy Majorana masses to vary in the TeV range. Thus mixing coefficients as 
large as advocated in Ref. {UeNjJ^N,) 2 — 10 _1 — 10~ 2 , could only arise in unnatural and 



fine tuned models 



221 ] . Note that, approximatively, the cross section goes like (s 



2 \2( 2 \2 



for real matrix elements. In this context, the coupling of HMN to gauge bosons and leptons 
is fixed to gUiNii where g is the SU{2) gauge coupling of the SM. Since the width of HMN 
grows with M^, at a certain value it will happen that Tjv > Mn, signaling a breakdown 
of perturbation theory. The perturbative limit on Mn is thereby estimated requiring < 
Mat/2, which gives an upper bound of ~ 3 TeV jl5. I^. 

In Fig.E^a) the cross section is plotted for masses up to this perturbative limit, using the 
maximally allowed value of the mixing. We see that for M^. = Mjv- = 3 TeV the signal does 
reach the level of 10 _1 , 10~ 2 fb respectively for the (rr) and the (fi/J,) signals at ^/s = 500 
GeV, which for an annual integrated luminosity of 100 fb -1 would correspond respectively 
to 10 and 1 event/year. At higher energies, O (TeV), one could get even larger event rates 
(30 and 3) respectively. The solid curve refers to e~e _ — ► r _ r _ : this is largest because the 
upper limits on the mixing are less stringent. One can also see the onset of the asymptotic 
regime at ^Js w 3 TeV. Fig. 0{b) shows that the cross section quickly decreases as lower 
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Major ana masses are considered. 

As even in the more optimistic cases event rates are quite modest it is important to 
check how the signal cross-section is affected by kinematic cuts on the angle of the outgoing 
leptons. The angular distributions turn out to be practically constant as shown in Fig. 0] 
They are forward-backward symmetric because both the t and u channel are present. The 
absence of a strong dependence on the polar angle is due to the fact that within the range 
of the parameters used here the contributing four point functions depend very mildly on the 
kinematic variables (u and t). This behaviour can be most easily understood using helicity 
amplitudes. The spinorial part common to all the diagrams is: 

[v(p 2 )P L u(p 1 )] [u(p 3 )Prv(p4)} = S(p 2 ,Pi)T(p 3 ,p 4 ), (10) 

that in the limit of massless external particles is a well defined helicity amplitude: e^L — > 
IlIl- In the center of mass frame this is a S-wave scattering with J z = 0, meaning that the 
scattered particles are emitted back to back but without a preferred direction relative to the 
collision axis (z). 

So this signal is characterized by practically flat angular distributions and as a result the 
total cross section is quite insensitive to angular cuts. Values of ot for different angular 
cuts are reported in Tabled With | cos#|<0.99 the change in ot is ~ 1% for all energies 
considered, while using | cos6*|<0.95 the total cross-section decreases by ~ 5%. Note that 
the reduction of the total cross section is measured almost precisely by the reduction of the 
phase space, meaning that the angular distribution is constant up to ~ 0.1%. Thus it can 
be concluded that the number of events will not be drastically affected for any reasonable 
choice of experimental cuts. 

In supersymmetric (SUSY) extensions of the see-saw framework (e.g. [3|) the natural 
mass scale of the singlet neutrino sector is at least of order O(10 12 ) GeV: in a unified 
scenario such a value improves the unification of gauge coupling constants [27] . Therefore 
HMN masses in the TeV range -although not ruled out experimentally- are disfavoured from a 
model building point of view. In the effective low energy SUSY see-saw framework, however, 
L is violated by the light SU{2)l doublet sneutrinos slQ, Q, S]- The mass states % 2 (£ 
denotes the generation) exhibit a mass-spi 
radiative contribution to neutrino masses 



itting Ami = — mt 2 that is constrained by its 
3fl| : Am; < 36(156)(m e e xp /leV) keV for a common 



scale of SUSY masses of 100 GeV, and m e * p the experimental limits on neutrino masses. 
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FIG. 4: Angular distribution in the polar angle of the outgoing lepton for different values of the 
center of mass energy, y/s in the case of e~e~ — * t~t~ with Mjv 4 = Mjv. = 3 TeV. The curves 
are not exactly constant, and using an appropriate scale they show small deviations from a stright 
line, remaining left-right symmetric. 
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The two different values refer to average and absolute upper limits when scanning over the 
SUSY parameter space. Then, in addition to the HMN mediated contribution discussed 
above, a diagram containing L- violating doublet sneutrinos and charginos is present, see 
Fig-EJe). The L- violating doublet sneutrino propagator and hence the resulting cross section 
is proportional to Am. Such a contribution has been considered in Ref. jsj for the case of eV 
scale sneutrinos. In the realistic case of (9(100) GeV scale sneutrinos, the exact differential 
cross section, in the notation of Fig. ^f ) for the momenta, is 



da 



dcosQ 
A 



/ a 



1287T \Xw ■ 

D 2 + j2D 2l . 

i=l 



A)o(M) + D 00 (s,u) 



uA(s,t) +tA(s,u) 



Here, the sum over the (maximally mixed) individual mass states in the L-violating sneu- 
trino propagators is included in the loop coefficients D 00 and A and the chargino is assumed 
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TABLE I: Total cross section (for two different angular cuts at some sample energies). The 
corresponding cuts on the transverse momentum of outgoing leptons are also shown. The numerical 
values for masses and mixing correspond to the solid line of Fig. 3(a) and Fig. 4. 

y/s = 300 GeV y/s = 500 GeV y/s = 800 GeV y/s = 1000 GeV 
<r(fb) a(fb) a(fb) a(fb) 
|cos0|<l 0.094 0.176 0.301 0.379 
| cos 6>|<0.99 0.093 0.174 0.297 0.374 
|cos0|<O.95 0.089 0.167 0.285 0.358 

to be a pure gaugino. In Fig. |S] the resulting total fi~ fi~ — > t~t~ cross section is plotted 
for a common SUSY mass (mi + 777,2) /2 = ml = m x = 100 GeV and for (maximal) values 
Am u = 30 GeV and Am T = 80 GeV (Amp must not exceed ml, otherwise the vacuum 

n n 

becomes unstable |2i|) allowed by the kinematic upper limits on neutrino masses [3Jj m^ < 
190 keV and m T < 18.2 MeV. Even for such unrealistically large neutrino masses and, cor- 
respondingly, unrealistic values for the sneutrino mass-splitting, the maximal cross section, 
at the threshold singularity for real x~X~ production, is of order (9(10 -3 ) fb and therefore 
too small to be observable even for a nominal integrated luminosity of 100 fb _1 /yr. The 
cross section with two colliding electrons is even smaller because mf p < 3 eV. Unlike the 
Majorana neutrino mediated contribution the SUSY cross section decreases for larger val- 
ues of ra[. Since the reaction e~e~ — ► i~t~ conserves total lepton number, processes like 
e~e~ — > e~ fi~ may arise due to interactions violating lepton flavour number but conserving 
the overall lepton number. This type of process will be discussed in the SUSY framework 
in a forthcoming paper . 

Concluding, we have calculated the cross section for the process e~e~ — > t~ t~ (t = fx, r) 
keeping the full dependence on the external momenta in the loop calculation and using the 
maximal value of effective light-heavy mixing angles allowed by experiments. We find that 
only for Majorana masses in the TeV range the reaction has a measurable cross section 
(above 10~ 2 fb) with better prospect for the r signal (the corresponding mixing coefficients 
being the .ess constrained), thereby arriving at somewhat .ess optimistic conc.usions than 
in Ref. [2|. We have also estimated the corresponding SUSY contribution arising from 
sneutrino mixing. This, although it is affected by an enhancement in the region of the 
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FIG. 5: Total cross section for the sneutrino mediated reaction fj,~ [i~ — * t~t~. The plot refers to 
the following choice of parameters: Am^ = 30 GeV and Am T = 80 GeV, m = m x = 100 GeV. 

threshold singularity, remains below the minimal observable value of 10 -2 fb even for the 
(unrealistic) maximal value of sneutrino mass splitting. 
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